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INFLUENCE  OF  THE  INTEGRATION  LIMITS  ON  THE  SHAPE  OF  PAIR 
CORRELATION  FUNCTIONS  OF  NON-CRYSTALLINE  MATERIALS 
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Influence  of  the  termination  effect  on  the  shape  of  pair  correlation  functions  and  values  of 
structural  parameters  of  liquid  and  amorphous  materials  has  been  investigated.  It  is  shown  that 
termination  of  the  diffraction  curve  at  Qm\n  > 0 rather  strongly  affects  the  region  around  the 
first  peak  of  g(r)  especially  its  left  flank.  On  the  other  hand,  presence  or  absence  of  a prepeak 
at  low  g-values  of  the  structure  factor  does  not  bring  any  difference  to  the  pair  correlation 
functions.  Decreasing  of  the  upper  limit  of  the  diffraction  vector  causes  shift,  decrease  and 
broadening  of  the  first  peak  of  the  pair  correlation  function.  Besides,  this  yields  spurious 
oscillations  in  the  vicinity  of  the  first  and  the  second  maximum  of  g(r)  of  the  liquid  and 
amorphous  materials.  Therefore  straightforward  interpretation  of  the  experimental  data  may 
be  misleading.  Some  procedures  to  solve  the  problem  are  considered. 
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1.  Introduction 

Studies  of  non-crystalline  structure  are  mainly  based  on  the  experimental  determination  of  the 
X-ray  or  neutron  diffraction  curve  as  a function  of  the  scattering  vector  Q(Q  = (47r/A,)sin$,  where  X is 
the  used  wavelength  and  0 is  half  the  scattering  angle).  Transformation  of  the  scattering  intensity 
7=7(0  to  the  pair  correlation  function  g(r)  and  (or)  to  the  radial  distribution  function  4nr2p(r ) yields 
such  structural  characteristics  as  the  nearest  neighbour  distance  J and  the  coordination  number  of  the 
first  coordination  sphere  as  well  as  corresponding  values  for  the  higher  coordination  spheres  which 
are  important  parameters  of  the  short-  and  medium-range  order. 

The  pair  distribution  function  g(r)  can  be  calculated  by  Fourier  transformation  of  the  S(Q) 
spectrum  using  the  standard  transformation  technique: 

g(r)  = = 1 + T-j—  Ie[S(0  - l]sin(e  ■ r)dQ  (1) 

Po  2?r  rp0  o 

Here  p(r)  and  p0  are  the  local  and  average  number  densities,  respectively,  and  S(Q)  is  the 
Faber-Ziman  total  structure  factor  [1]  obtained  from  the  coherent  scattering  intensity: 


5(0  = 


::*(0-{{/2(0)-(/(0)2} 

</(0)2 


(2) 


with 


(/2(0)  = X erf®),  (AQj)  = X Cjfj  (0, 


(3) 


where  cj  is  the  mole  fraction  and  f/Q)  the  total  atomic  scattering  factor  of  the  /-th  component  in  the 
alloy. 
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The  position  of  the  first  maximum  of  g(r)  is  the  mean  nearest  neighbour  distance  /.  The 
coordination  number 


N = \Am'2  p{r)dr 

n 

is  the  mean  number  of  atoms  within  the  coordination  sphere  limited  by  /•/  and  The  first 
coordination  number  can  be  estimated  by  several  methods.  Most  commonly  used  ones  are:  i) 
symmetrical  (the  coordination  shell  is  considered  as  symmetrical  about  a first  maximum  in  the  g(r) 
curve,  i.e.  the  right-hand  side  of  the  peak  is  made  symmetrical  with  that  on  the  left)  and  ii)  integration 
from  the  left-hand  edge  of  the  first  peak  to  the  first  minimum  in  the  radial  distribution  function. 

The  transform  in  Eq.(l)  requires  data  from  Q = 0 to  Q = <x>.  In  principle,  the  scattering 
intensity  can  be  measured  down  to  Q - 0.  However,  experimental  limitations  and  in  some  cases 
additional  small  angle  scattering  lead  to  a lower  limit  Qm\n  > 0 for  each  diffraction  experiment.  On  the 
other  hand,  the  diffraction  vector  is  terminated  by  a certain  value  Qmm  which  is  determined  by  the 
wavelength  used:  Qmm  = 4ji/A,.  As  a result,  the  integrand  in  Eq.(l)  is  limited  by  Qmm  and  0„,ax,  that  is: 


1 £?max  r 

g"'(r)  = 1+— f-  \Q[S{.Q)-\\MQ-r)dQ 

£71  0min 


(5) 


Due  to  such  limitation  of  the  diffraction  vector,  the  resulting  pair  distribution  function  gexp(r)  differs 
from  the  ideal  g(r).  So-called  termination  or  truncation  effect  yields  spurious  details  on  the  pair 
correlation  functions.  This  deforms  the  radial  distribution  functions,  changes  the  position  and  profile 
of  the  maxima  and  minima  and  introduces  errors  in  the  calculated  coordination  numbers  [2  - 5]. 

Jn  this  paper  we  examine  influence  of  the  termination  effect  on  the  shape  of  the  pair 
correlation  functions  of  the  non-crystalline  materials  with  low-coordinated  (open)  structure.  Liquid  Te 
and  amorphous  Ge-telluride  have  been  chosen  for  the  investigation. 


2.  Experimental 

Neutron  scattering  experiment  on  liquid  Te  has  been  carried  out  on  the  double-axis 
spectrometer  of  the  research  reactor  at  Rossendorf.  The  neutron  wavelength  was  A.  = 0.91  A and  the 
accessible  0-values  range  — from  1.0  A’1  to  9.35  A'1.  The  sample  was  contained  in  a sealed  quartz 
tube  (internal  diameter  15  mm,  wall  thickness  1 mm).  The  overall  uncertainty  in  the  sample 
temperature  was  smaller  than  5 K.  The  data  correction  procedure  has  been  described  in  [6]. 

The  X-ray  diffraction  experiments  were  performed  with  a 6-6  diffractometer  using  MoKa~ 
radiation  which  has  been  selected  by  a focusing  graphite  monochromator  placed  in  the  diffracted 
beam.  The  X-ray  intensity  scattered  at  the  free  surface  was  measured  by  a scintillation  detector  with  a 
pulse  height  analyser.  The  liquid  Te  temperature  was  constant  within  an  error  of  about  1 K.  The 
scattering  intensity  has  been  measured  from  1 .0  to  1 1 .85  A’1  of  Q in  the  case  of  liquid  Te. 

The  amorphous  Ge25Te75  alloy  was  prepared  by  the  melt-spinning  technique.  The  X-ray 
intensity  diffracted  by  the  amorphous  sample  at  room  temperature  was  counted  over  the  range  of  Q 
between  0.5  and  12.5  A'1.  The  correction  and  normalization  of  the  experimental  data  were  done  as 
described  in  [7  -9]. 
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3.  Results  and  discussion 

3.1.  Liquid  Tellurium 

The  experimental  structure  factors  of  liquid  Te  at  550  °C  (neutron  diffraction)  and  500  °C  (X- 
ray  diffraction)  are  shown  in  Fig.  1.  The  differences  between  the  results  of  X-ray  and  neutron 
scattering  are  mainly  caused  by  the  rather  sensitive  absorption  correction  of  the  neutron  experiment. 
In  both  cases  scattered  intensity  has  been  measured  down  to  Qm = 1.0  A'1.  For  smaller  g- values  we 
have  employed  a parabolic  extrapolation  of  S(Q): 

S(Q)  = S(0)  + AQ2.  (6) 

The  value  of  S( 0)  = 0.062  has  been  taken  from  Ref.  [10]  and  the  coefficients  was  determined 
by  fitting  the  extrapolation  function  to  the  experimental  S(Q)  at  Q = Qmm- 
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Neutron  diffraction,  550  °C 
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Fig.  1 . Structure  factors  of  the  liquid  Te:  circles  - experimental  results;  lines  - extrapolation. 

Fig.  2 shows  the  pair  correlation  and  radial  distribution  functions  calculated  with  Eq.  (5)  using 
lower  integration  limits  of  Qm in  = 0 A'1  and  Qmm  = 1 .0  A'1.  The  upper  integration  limit  was  Qmax  = 9.35 
A'1  for  neutrons  and  Qmax  = 1 1.85  A-1  for  X-rays.  It  is  very  good  seen  that  the  region  around  the  first 
coordination  shell,  especially  its  left  flank,  is  rather  strongly  affected.  Termination  of  the  structure 
factor  at  Qmin  > 0 substantially  increases  the  pair  correlation  and  radial  distribution  function  in  the 
region  of  the  first  maximum  that  yields  overestimated  coordination  number. 


4nr  p(r) 
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Fig.  3.  Pair  correlation  function  of  liquid  Te  calculated  with  different  integration  limits  Qimx . 


Fig.  4.  Residual  pair  correlation  function  A g(r)  calculated  with  S(Q)  from  gmax  to  oo. 
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The  next  step  was  to  examine  the  effect  of  the  structure  factor  termination  at  high  0-values. 
For  this  purpose  we  have  extrapolated  the  experimental  structure  factor  for  Q > Qmm  using  the 
function 


S(Q)  = 1 + c,  cos(e,0  - c3)exp(-c40)/0  (7) 

which  enables  one  to  obtain  an  analytical  expression  for  g(r)  up  to  infinity.  The  parameters  C1...C4 
have  been  determined  by  a least  squares  fit  of  the  high  0-region  of  the  experimental  S(Q).  The 
extrapolated  parts  of  the  liquid  Te  structure  factors  are  given  in  Fig.  1. 

In  Fig.  3 the  region  of  the  first  and  the  second  maximum  of  the  pair  correlation  function  of 
liquid  Te  at  500  °C  calculated  with  different  upper  0-limit  is  shown.  It  is  seen  that  with  increasing 
0max  the  height  of  the  first  peak  increases  while  its  width  decreases.  Besides,  the  position  of  the  first 
maximum  is  shifted  to  the  smaller  r-values  and  “false”  oscillations  disappear. 

Fig.  4 illustrates  the  termination  error  in  g(r)  by  limiting  the  structure  factor  to  0max  = 9,  13, 
17  and  21  A'1.  It  is  seen  that  the  high-angle  contribution  to  the  pair  correlation  function 

Ag(r)  = 1+—I—  ]Q[S(Q)-\]sm(Qr)dQ  <8) 

27r  fp0  Qnux 

is  the  most  essential  in  the  vicinity  of  the  first  and  the  second  peak. 

The  pair  correlation  functions  of  the  liquid  Te  calculated  with  different  lower  and  upper 
integration  limits  in  Eq.  (5)  are  presented  in  Fig.  5.  It  is  seen  that  both  limitations  significantly  affect 
the  resulting  curve,  especially  in  the  range  of  the  first  and  the  second  coordination  spheres.  Therefore 
a straightforward  interpretation  of  the  “experimental”  pair  correlation  functions  may  be  misleading. 


2 3 4 5 6 7 8 

r,  A 

Fig.  5.  Pair  correlation  functions  of  the  liquid  Te  calculated  with  different  integration  limits 

0imn  and  0 max- 


In  Fig.  6 we  show  some  structural  parameters  of  liquid  Te  in  dependence  on  the  upper 
integration  limit  0maX.  The  nearest  neighbour  distance  r and  coordination  number  N1  continuously 
decrease  with  increasing  0max  up  to  19-20  A'1  while  the  value  of  g(7)  increases.  Beginning  with 
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<2max  ~ 20  A'1  all  parameters  remain  practically  unchanged.  This  is  in  good  correspondence  with  the 
result  shown  in  Fig.  4.  Indeed,  the  residual  pair  correlation  function  A g(r)  becomes  very  small  when 
gmax  achieves  ~ 21  A'1  or  more. 


Fig.  6.  Nearest  neighbour  distance  /,  coordination  number  N'mm,  Nrsym  and  g(r')  as  a function 

Of&na*(0min  = O). 


3.2.  Amorphous  Ge2sTe75 

The  experimental  structure  factor  of  the  amorphous  Ge25Te75  alloy  measured  at  room 
temperature  is  shown  in  Fig.  7.  It  is  seen  that  the  modulations  of  the  S(Q)  of  the  amorphous  Ge- 
telluride  are  much  more  pronounced  than  those  of  liquid  objects  (e.  g.  liquid  Te).  Besides  the  shape  of 
the  maxima  and  minima  in  the  S(Q)  of  the  a-Ge25Te75  alloy  are  much  more  complex  than  that  in  the 
hard-sphere  structure  factor.  This  means  that  one  can  not  analytically  extrapolate  the  experimental 
structure  factor  in  the  same  or  similar  manner  as  we  did  it  for  the  liquid  Te.  In  such  cases  it  is  better  to 
use  experimental  data  up  to  the  limited  value  of  the  diffraction  vector  Qmwi  than  to  introduce 
additional  errors  employing  a theoretical  function  for  Q > Qmm. 


Fig.  7.  Experimental  structure  factor  of  the  amorphous  Ge25Te75  alloy. 

As  we  have  seen  above,  the  termination  of  the  structure  factor  at  Q - Qmax  yields  ghost 
oscillations  in  the  vicinity  of  the  first  peak  of  the  pair  correlation  function.  It  is  very  interesting  and 
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important  to  know  whether  these  false  modulations  depend  on  the  value  of  Qmm  with  regard  to  the 
shape  of  the  structure  factor.  In  order  to  test  it,  we  have  calculated  the  residual  pair  correlation 
function  using  Eq.  (8)  with  different  limit  Qmax.  As  a structure  factor  we  employed  the  analytical 
function  given  by  Eq.  (7)  with  parameters  ct  = 0.9239,  c2  = 2.8576,  c3  =20.2316,  c4  =0.2012.  Figs.  8 
and  9 show  the  high-angle  contribution  to  the  pair  correlation  function  A g(r)  in  dependence  on  gmax. 
The  oscillations  of  the  residual  pair  correlation  function  are  the  greatest  when  the  structure  factor  is 
terminated  at  the  position  where  the  S(Q)  equals  unity  and  shows  maximum  slope  (e.g.  at  13.1  A'1  and 
14.2  A'1).  On  the  other  hand,  oscillations  of  the  Ag(r)  are  minimal  when  Qmm  falls  at  the  position  of  a 
minimum  or  a maximum  (e.g.  at  13.7  A'1)  at  the  end  of  the  structure  factor  where  the  slope  of  S(Q)  is 
zero. 


Fig.  8.  Residual  pair  con-elation  function  A g(r)  calculated  with  S(Q)  from  Qm„  to  oo. 


Fig.  9.  Ag(71|r,0  as  a function  of  gmax  in  Eq.(8).  S(Q)  - part  of  the  structure  factor 
of  liquid  Te  at  500  °C  (as  in  Fig.  1). 


Taking  into  account  the  result  mentioned  above  we  have  truncated  the  experimental  structure 
factor  of  the  a-Ge25Te75  alloy  at  Qmilx  of  1 1 .2  and  13.5  A'1  that  respects  to  the  position  of  the  minimum 
at  the  end  of  curve.  Fig.  10  illustrates  that  decrease  of  the  upper  integration  limit  in  Eq.  (5)  causes 
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decrease,  shift  and  broadening  of  the  first  peak  and  yields  spurious  oscillations  in  the  pair  correlation 
functions  of  the  amorphous  material  in  the  same  manner  as  in  the  case  of  liquids. 


Fig.  10.  Pair  correlation  (a)  and  radial  distribution  (b)  functions  of  amorphous 
Ge25Te75  alloy  calculated  with  different  integration  limits  gmax. 


a) 


b) 


Fig.  1 1.  Pair  correlation  (a)  and  radial  distribution  (b)  functions  of  amorphous  Ge25Te75  alloy. 


In  the  experimental  structure  factor  of  the  a-Ge25Te75  alloy  (Fig.  7)  one  can  see  a small 
prepeak  at  Q = 1.1  A'1.  It  is  generally  assumed  that  i)  this  feature  of  the  diffraction  pattern  is 
associated  with  medium-range  order  ( MRO ) and  (or)  ii)  prepeak  reflects  strong  bonding  of  unlike 
atoms.  As  regards  the  amorphous  Ge25Te75  alloy  we  believe  that  the  prepeak  in  the  structure  factor 
arises  due  to  connection  of  Ge  and  Te  atoms.  However  in  the  present  paper  we  would  like  to  examine 
whether  and  how  the  shape  of  the  pair  correlation  function  depends  on  the  prepeak.  Therefore  we 
calculated  pair  correlation  function  using  experimental  structure  factor  i)  from  Qmin  = 0.5  A'1  to  Qmm  = 
13.5  A'1;  and  ii)  from  Qm j„  = 1.5  A'1  to  Qmm  = 13.5  A'1.  In  both  cases  S(Q)  was  extrapolated  from 
SfOmm)  to  S{ 0)  as  described  above  for  the  liquid  Te.  As  one  can  see  from  Fig.  1 1 there  is  practically 
no  difference  between  pair  correlation  functions  calculated  with  structure  factor  containing  prepeak 
and  without  that. 
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4.  Conclusions 

Limitation  of  the  experimental  scattering  intensity  at  certain  values  of  the  diffraction  vector 
Qm in  and  Qmax  that  differ  from  zero  and  infinity,  respectively,  yields  spurious  details  on  the  pair 
correlation  functions  and  introduces  errors  in  the  calculated  structural  parameters  of  the  liquid  and 
amorphous  materials. 

Termination  of  the  diffraction  curve  at  Qm\n  > 0 significantly  increases  the  pair  correlation  and 
radial  distribution  functions  in  the  region  of  the  first  maximum  that  yields  overestimated  coordination 
number. 

Termination  of  the  structure  factor  at  high  (Lvalues  of  the  diffraction  vector  causes  shift, 
decreasing  and  broadening  of  the  first  peak  of  the  pair  correlation  functions.  In  addition,  this  results  in 
spurious  oscillations  in  the  vicinity  of  the  first  and  the  second  maximum  of  g(r). 

Termination  effect  in  the  pair  correlation  functions  can  be  eliminated  or  significantly  reduced 
by  extrapolation  of  the  experimental  scattering  intensity  or  structure  factor  to  zero  and  to  infinity 
using  analytical  functions.  In  cases  when  this  cannot  be  done  the  termination  at  high  g-values  should 
be  made  at  the  position  of  minimum  or  maximum  of  the  structure  factor. 
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